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ABSTRACT 

We present maps of 7.78 square degrees of the Lupus molecular cloud complex at 24, 70, and 160/im. 
They were made with the Spitzer Space Telescope's Multiband Imaging Photometer for Spitzer (MIPS) 
instrument as part of the Spitzer Legacy Program, "From Molecular Cores to Planet-Forming Disks" 
(c2d). The maps cover three separate regions in Lupus, denoted I, HI, and IV. We discuss the c2d 
pipeline and how our data processing differs from it. We compare source counts in the three regions 
with two other data sets and predicted star counts from the Wainscoat model. This comparison shows 
the contribution from background galaxies in Lupus I. We also create two color magnitude diagrams 
using the 2MASS and MIPS data. From these results, we can identify background galaxies and 
distinguish them from probable young stellar objects. The sources in our catalogs are classified based 
on their spectral energy distribution (SED) from 2MASS and Spitzer wavelengths to create a sample 
of young stellar object candidates. From 2MASS data, we create extinction maps for each region and 
note a strong corresponence between the extinction and the 160 /im emission. The masses we derived 
in each Lupus cloud from our extinction maps are compared to masses estimated from ^^CO and C^^O 
and found to be similar to our extinction masses in some regions, but significantly different in others. 
Finally, based on our color-magnitude diagrams, we selected 12 of our reddest candidate young stellar 
objects for individual discussion. Five of the 12 appear to be newly-discovered YSOs. 
Subject headings: infrared: stars — ISM: clouds — stars: formation 



1. INTRODUCTION 

The Spitzer Legacy Science Program, "F rom Molecu- 
lar C ores to Planet-Forming Disks" (c2d) () Evans et al.l 
l2003l ). mapped five nearby molecular clouds us- 
ing the Spitzer S pace Telescope (hereafter Spitzer) 
(jWerner et al.ll2004[ ). This paper, fourth in a series of pa- 
pers, presents the basic results from the Multiband Im ag- 
ing Photometer for Spitzer (MIPS) (iRieke et al.ll2004f l for 
the Lupus Molecular cloud complex. 
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The Lupus molecular cloud complex covers a broad 
range of Galactic longitude and latitude and consists of 
a number of distinct regions (|Krautteij|199l[ ). The cloud 
is located near the Scorpius-Centaurus OB association. 
The large spatial extent of Lupus on the sky (~ 20°) 
means that different regions in the cloud may be at dif- 
ferent distances. Hipparcos measurements to nearby OB 
star sub-groups within Scorpius-Centaurus g ive average 
distances of ~ 140 pc (jde Zeeuw et al.lfl999t ). This dis- 
tance is in ag reement with various reddening studies (e.g. 
iFrancd 1990). Distance estimates derived from Hippar- 
cos parallaxes to individual stars within the Lupus clouds 
can vary by more than a factor of two, suggesting that 
depth effects are important in the Lupus clo ud complex 
(jWichmann et all 119981 : iBertout et all 119990 . For this 
paper we are adopting distances of 150 ±20 pc for Lupus 
I and IV and 200 ± 20 pc for Lupus HI. See IComeroiil 
(|2007l and references therein) for a full discussion on the 
difficulties in obtaining accurate distance measurements 
to Lupus. 

Compared to other c2d cl ouds. Lu pus is a relatively 
quiescent star-forming cloud. IComer on (2007) has com- 
piled a table of all the known classical T-Tauri Stars 
(cTTSs) in different Lupus regions. Within the area we 
observed with bot h the InfraRed Array Camera (IRAC) 
(|Fazio et al.l l2004[ ) and MIPS instruments, there are 5 
CTTSs in Lupus I, 32 in Lupus III, and 3 in Lupus 
IV. In addition to the known CTTSs, IComeronI ()2007D 
lists 19 new suspected very lo w-mass members of Lupus 
HI with likely Ha emission (|L6pez Martf et al.l [2005), 
and 4 recently discovered M-type members in Lupus III 
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TABLE 1 
Summary of Observations 
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(jCoineron et alJl2003f ). Over 100 weak-li ne T-Tauri stars 
(WT TSs) were discovered with ROSAT ("Kr autter efan 
Il997| ). extending over a much larger region than the 
dark clouds. Based on v elocity dispersions and positions, 
IWichmann et all (|1997f) speculated that these WTTSs 
may have formed earlier and since moved away from the 
dark clouds. 

We summarize our observations in § 2 and our data 
processing pipeline in § 3. A full description of the data 
processing used on all c2d data is av ailable in our delivery 
documentation (jEvans et al.ll2007b . For this paper, we 
made 'MIPS high-reliability' cuts to the c2d catalogs, as 
described in § 4. Our results are summarized in §[5l 

Our goals in this paper are to present the c2d Lupus 
MIPS data and identify the likely young stellar mem- 
bers of the Lupus cloud. We plot the differential source 
counts at 24 /im for our cloud regions and compare these 
to model predictions and observed off-cloud fields (§[5Tl|). 
We use the 2MASS and MIPS data to create two color- 
magnitude diagrams (§[5?2]) which are useful in distin- 
guishing background galaxies from young stellar objects. 
We create a sample of 'Young Stellar Object (YSO) can- 
didates' based on each source's spectral energy distribu- 
tion in the 2MASS and Spitzer wavebands and discuss 
the significance of our sample in relation to previously 
know YSOs fij l5.3p . Because YSOs are typically found 
in the densest regions, we create extinction maps from 
2MASS, estimate the mass in each cloud region, and 
compare the extinction maps to the 160 fim emission 
(i) l5.4|) . Finally, from our 'YSO candidate' sample and 
the color-magnitude diagrams, we select 12 likely young 
stellar objects. These sources are discussed individually 
in §[5111 

2. OBSERVATIONS 

Our observations of Lupus cover three separate on- 
cloud regions and two off-cloud regions. We designate 
the on-cloud regions as I , III, and IV following the con- 
vention from lCambresvl (|1999D . The total mapped area 
on-cloud in two epochs is approximately 7.78 square de- 
grees. In addition to the on-cloud maps, we observed 
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Fig. 1. — Map showing the observed regions in the Lupus cloud. 
These regions are designated: Lupus I, Lupus III, Lupus IV, and 
two off-cloud regions (OCl and OC2). The outlines show the 
mapped areas at 24 fim; the 70 and 160 fim coverage maps are 
very similar. The grayscale image is the IRAS 25 fim emission. 
The saturated region at the lower left is the Galactic plane. The 
directions of increasing Galactic longitude and latitude are denoted 
by arrows in the lower right corner. 

two small off-cloud regions for reference, designated off- 
cloud 1 (OCl) and off-cloud 2 (0C2). In Figure [H the 
observed 24 fim areas for each region are overlaid on the 
25 /im IRAS map. The observed areas at 70 and 160 iim 
are very similar in size, shape, and position to those at 
24/im. Tablc[T]lists the Astronomical Observing Request 
(AOR) numbers, dates observed, and Spitzer program ID 
for all our observations. 

We used the optical extinction maps of ICambresvl 
()1999( l to choose the areas mapped. Lupus I covers ar- 
eas with extinctions Ay > 3, while Lupus III and IV 
map extinctions Ay > 2. The two off-cloud regions had 
Ay < 0.5 in the optical extinction maps. Lupus II was 
not observed because there was very little area in the 
extinction map with Ay > 2. Due to time constraints, 
we did not observe Lupus V or VI. It is worth noting, 
however, that there is no obvious evid ence of recent o r 
current star formation in these regions (jComero 
Information on the observed regions, including mapped 
area, distance, and approximate Galactic coordinates are 
hsted in Table El 

All data were taken in MIPS fast-scan mode with a 
240" offset between scan legs. A second epoch of obser- 
vation was taken later to help identify asteroids. The 
time between the two epochs of observation ranged from 
3.5 hours in Lupus I to almost 7 hours in Lupus III. The 
second epoch is offset, compared to the first, by 125" in 
the cross-scan direction. Furthermore, to fill in gaps at 
160 fim, the second epoch is offset by 80" along the scan 
direction. Only half of the 70 ^m array returns usable 
data. To minimize this effect, the AORs were planned 
such that the second epoch would fill in the holes in the 
70 fim maps. The net result is that the 24 and 70 fim 
maps have approximately the same coverage, but the to- 
tal integration time per pixel at 70 fim is 15 seconds, 
compared to 30 seconds at 24 /im. The integration time 
per pixel is roughly 3 seconds at 160 fim. 

3. DATA REDUCTION 
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TABLE 2 

Basic Information on the Observed Regions 



Region 


/ 

(deg.) 


b 

(deg.) 


dist. 

(pc) 


size 
{dcg2) 


median Ay 
(mag) 


Mass (Ay > 3) 
(Mq) 


Lupus I 


338 


17 


150 


3.82 


1.8 


440 


Lupus III 


340 


9 


200 


2.88 


2.4 


690 


Lupus IV 


336 


8 


150 


1.08 


2.2 


250 


Lupus OCl 


344 


12 




0.31 


0.7 




Lupus OC2 


343 


4 




0.31 


2.1 





Note. — The approximate Galactic latitude and longitude, distance, 
mapped area, median Ay, and mass are listed for the thr ee on -clo ud regions 
and th e two off-cloud regions. The masses are derived in ^ 15.41 See IComerorJ 
120071) for a discussion of the adopted distances to each region. The map size 
listed is the size of the 24 /xm map covered in both epochs. The regions observed 
at 70 and 160 fim are of similar size. 



The c2d team has made regular dehveries of mosaics 
and bandmerged catalogs to the Spitzer Science Center 
(SSC) for c2d regions. The delivery documentation and 
data used in this paper can be found on the SSC web- 
site under delivery 4, the final data delivery of the c2d 
team^ . Detailed discussion of the data reduction method 
applied to all c 2d products c an be found in the delivery 
documentation (jEvans et al . 2007). In this paper we will 
give only a brief outline of the standard processing pro- 
cedure. We will discuss the steps of this procedure in the 
order they are performed. For this paper, as in delivery 
4, we processed all the data with the S13 SSC pipeline. 

We made two changes to the standard c2d processing 
procedure. First, for display purposes only, we made 
several cosmetic improvements to the 70 fim and 160 fim 
images. These enhanced images were created from the 
older Sll SSC products. The specific improvements will 
be discussed in ^ \3.1\ Source extraction and photome- 
try at these wavelengths is performed with the standard 
SSC filtered and unfiltered data products using the c2d 
pipeline. Our second change compared to the standard 
c2d processing procedure is that at 160 /im we band- 
merged the c2d extractions with the c2d catalog. 

This paper is one of a series of papers to present the c2d 
observations of clouds in a uniform manner. The first pa- 
per in__tliis_series presented MIPS results for Chamaeleon 
II (I Young et al l lBol . That paper used an earlier ver- 
sion of the SSC pipehne, S9.5. Other papers that present 
MIPS results fro m the c2d c louds use t he S13 products 
as we have he re (iRebull et al . 2007; Har vev et ani2007l : 
iPadgett et aI]|2007D . Despite some differences in data 
processing, many of the figures are standardized to facil- 
itate comparisons between c2d clouds. 

3.1. Image Processing 

Starting from the Basic Calibrated Datasets (BCDs), 
the c2d team further processed the 24 fim images to re- 
move some instrumental effects. We corrected a "jailbar" 
response pattern that is caused by bright sources or cos- 
mic rays by applying an additive correction to the low 
readouts to bring them up to the same level as the high 
readouts. Additionally, the first four frames in a map 
were scaled to the median of following frames, and each 
AOR was self-flattened to eliminate a small instrumental 
brightness gradient along the column direction and other 
residual image artifacts. 

^ http:/ /ssc. spitzer. caltech.edu/legacy/c2dhistory.html 



The SSC provided two sets of BCD products for the 
70 and 160 fim data, unfiltered (normal processing) and 
filtered (a time median filter is applied) . We have further 
processed the SSC's unfiltered 70 /im data to improve 
the visual appearance of the mosaic image. The first 
correction we made was to subtract out the latent signal 
from stimulator (stim) flashes. The signal after the latent 
had faded was used as a reference to subtract out the 
stim flash latents. Further, the SSC unfiltered products 
contained a striped pattern roughly along the column 
direction, which was caused by response variations. We 
self-flattened the data with a median obtained from the 
frames in each scan leg to reduce this effect. 

We have made visually enhanced 160 /im mosaics by 
excluding the stim flash BCDs and the following frames 
that contained a latent signal. The mosaic was then cre- 
ated with 8" pixels, half the natural pixel scale (twice the 
resolution). Lastly, the image was smoothed with a me- 
dian filter, where each pixel in the image was replaced by 
the median of the values within a 7 x 7 pixel box. Pixels 
with a value of NaN were treated as missing data during 
this filtering. 

Figures [2] - [4] show the 24,70, and 160 /im images in 
separate panels for regions I, III, and IV, respectively. 
We have also combined these three wavelengths to make 
a color image where red represents the 160 /im, green the 
70 /xm, and blue the 24 fim emission. The color images 
for Lupus I, III, and IV are shown in Figures [5] - [T) 

3.2. Mosaicking and Source Extraction 

We created mosa ic images at 24 fim using the SSC's 
MOPEX program (jMakovoz fc Marleaul [20051 ) with the 
outlier rejection and position refinement modules. Three 
sets of mosaics were created, one for each epoch sepa- 
rately and one combined epochs mosaic. We used the im- 
proved BCDs described in § l3.1l for these mosaics. Outlier 
rejection excluded bad pixels caused by cosmic rays or 
image artifacts. The initial identification of point sources 
was performed on each mosaic using a modified version 
of the DOPHOT program (Schcchter et al. 1993). This 
program works by fitting point spread functions, so as 
sources were found, they were subtracted from the im- 
age. This procedure was iterated until a user defined flux 
limit was reached. We calculated the final position and 
photometry for each source using all BCD images that 
contributed to that source. 

The 24 fim source fiuxes were calibrated by comput- 
ing the ratio between the aperture flux and model point 
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Lupus I 




RA (JSOOO) 

Fig. 2. — The 24, 70, and 160 images of Lupus I. We have applied several cosmetic corrections to the 70 and 160 fim images as 
described in ii l3.1l 



Lupus in 




RA (J3000) 

Fig. 3. — The 24, 70, and 160 /im images of Lupus III. We have applied several cosmetic corrections to the 70 and 160 fj,m images as 
described in S IS. II 
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Lupus IV 
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Fig. 4. — The 24, 70, and 160 fira images of Lupus IV. We have apphed several cosmetic corrections to the 70 and 160 /im images as 
described in 5 13. II 
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Fig. 5. — A pseudo-color image for Lupus I showing the 24, 
70, and 160 /im emission in the blue, green, and red channels, 
respectively. We have applied several cosmetic corrections to the 
70 and 160 fim images as described in ^ 13.11 
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Fig. 6. — A pseudo-color image for Lupus III showing the 24, 
70, and 160 iim emission in the blue, green, and red channels, 
respectively. We have applied several cosm etic corrections to the 
70 and 160 fj,m images as described in ^ 13.11 

spread function (PSF) flux for all sources with a reliable 
aperture flux that could be fitted by the point source 
profile. We assumed the average aperture to PSF flux 
ratio held for all sources, and a multiplicative correction 
was applied. 

The 70 /im point source extractions were performed 
on mosaic images made using MOPEX on the SSC's fil- 
tered and unfiltered BCDs, not the enhanced mosaics de- 
scribed in ^ 13.11 We used APEX, the SSC's point source 
extraction software, on both the filtered and unfiltered 
mosaic images to obtain point sources. We used the point 
response function (PRE) photometry from the filtered 
mosaic except for sources brighter than ^ 2 Jy, when 
the SSC's filtering over corrects in the wings of the PRE. 
In such cases, we substituted PRE photometry from the 
unfiltered mosaic. 
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Fig. 7. — A pseudo-color image for Lupus IV showing the 24, 
70, and 160 /im emission in the blue, green, and red channels, 
respectively. We have applied several cosmetic corrections to the 
70 and 160 fim images as described in ^ 13.11 

At 160 fj,m we performed source detection and pho- 
tometry using the MOPEX package on the SSC's un- 
filtered BCD mosaics. Again, we did not use the en- 
hanced mosaics created in § 13.11 Eluxes were determined 
by point-source fitting. Sources that were close multiples, 
extended, or had S/N < 5 are not hsted here. Eor satu- 
rated sources, we assumed the object was a point source 
and fit the wings of the PSE to find the flux. Sources 
which were detected but deemed unreliable by visual in- 
spection were dropped. Because of this culling, the 160 
/xm source list is believed to be reliable but not complete, 
and no completeness limit is given. 

The c2d team estimates the absolute flux uncertainty 
is 15% at 24/zm and 20% for 70 and 160/im. We estimate 
the relative flux uncertainty to be approximately 5% plus 
the statistical uncertainty for each source. 

3.3. Bandmerging 

The final stage of data processing prior to science anal- 
ysis of c2d images involves the creation of bandmerged 
catalogs. Again, refer to the c2d delivery document 
(jEvans et al.ll2007l ) for details. Even though this paper 
focuses on MIPS results, IRAC bandmerging is also de- 
scribed since these wavelengths are used for source clas- 
sification (§ 13. 4p and when discussing individual sources 

3.3.1. IRAC - MIPS 24 Bandmerging 

Eor each band, the three source extraction lists 
(epochl, epoch2, and combined epochs) were checked for 
"self-matches" within an epoch; two sources extracted 
within one epoch, but with positional matches of < 2'.'0 
were considered to be the same source. The fluxes of the 
detections were summed, and the position of the source 
was calculated as the weighted mean. We then merged 
the three lists together to cross-identify sources with po- 
sitional matches of 2'.'0 or less. We visually inspected 
the epoch-merged source lists for each band to remove 
diffraction spikes, column pull-down, latent images, and 
other image artifacts that were misidentified as sources. 
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The epoch-merged source lists for each band were then 
merged as follows: First, we merged the four IRAC bands 
together, one-at-a-time. We started by merging IRACl 
(3.6 /im) and IRAC2 (4.5 ^m), then merged this product 
with IRAC3 (5.8 /im), and finally combined with IRAC4 
(8 /im). At each step, we combined detections at two 
wavelengths into a single source if the difference in cen- 
tral position was < 2'.'0. Then, this bandmerged IRAC 
catalog was merged with the 24 /xm MIPSl band using 
a larger distance, 4'.'0. The larger distance was used be- 
cause of the larger PSF at 24 nm compared to IRAC. 
Lastly, the IRAC-l-MIPSl catalog was compared with 
the 2MASS catalog using a position matching criterion 
of 2'.'0. 

In all steps of bandmerging, IRAC, MIPSl, and 
2MASS, we merged together the closest position match 
between 2 wavelengths. Any other detections within the 
specified radius were preserved in the final catalog, but 
not merged. 

3.3.2. MIPS 70 nm and 160 /im Bandmerging 

Even though we have two epochs of observation at 
70 /im, they do not have much overlap between them. 
Thus, we only performed source extractions on the com- 
bined epochs dataset. First, we performed a 'self- 
merging' on the combined epochs dataset using a match- 
ing radius of 4'.'0. Secondly, each 70 /xm point source 
detection was matched to all shorter wavelength catalog 
sources within a S'.'O radius. We then visually inspected 
each merged 70 /im detection to determine the 'best' 
source match between shorter wavelengths (1.25 — 24/im) 
and 70 /im. This 'best' source match was the one with 
the most consistent Spectral Energy Distribution (SED) 
across all detected wavelengths. In a few cases where it 
was not clear which SED was 'better', we chose the closer 
match. 

Shorter wavelength sources that were within S'.'O of a 
70 /im detection but not matched with the 70 /im, were 
assigned a special flag to denote their status. Thus, in 
our final catalog, a 70 /im detection is listed once; we did 
not try to split up the flux to different shorter wavelength 
sources. 

There are so few 160 /im point sources that they were 
bandmerged individually by eye for those selected sources 
described in §[531 

3.4. Source Classification 

After bandmerging, we classified sources based on their 
SEDs at wavelengths from 1.25 — 24 /im. Even though 
this paper concentrates on the MIPS results for Lupus, 
we utilized the IRAC and 2MASS wavelengths for source 
classification. The two most important categories for the 
work discussed here are 'star' and 'Young Stellar Object 
candidate (YSOc)'. Objects classified as 'star' have mea- 
surements in 3 or more bands and the SED can be fitted 
by a reddened stellar photosphere (including the 2MASS 
data if available). In addition, if a source had an SED 
that could be fitted by a reddened photosphere except in 
one band, we also classified it as 'star', assuming that it 
has a prominent spectral line feature that changed the 
shape of the SED. A future pap er will explore our source 
classification method in depth f Lai et alllin prepj ). 

It is important to include the effect of dust extinction 
when identifying stars, because star light is often dimmed 



by significant dust in the molecular clouds we observed. 
The observed flux of a reddened star, FotsW, can be 
described with the following equation: 

l0g(F„fc,(A)/Fwei(A)) = log(X) - 0.4 X Ce.t(A) X Ay 

(1) 

where FmodeiW is the stellar photosphere model, K is 
the scaling factor of the model for a particular star, and 
CextW = ^a/^f is the ratio of extinction at wavelength 
A to visual extinction from the dust extinction law. K 
and Ay can be derived from the linear fit of this equa- 
tion by adopting stellar photosphere and dust extinction 
models. The stellar photosphere models for i^Tg-MIPSl 
bands are based on the Kurucz-Lejeune models and come 
from the SSC's "Star-Pet" tooP. For the 2MASS bands, 
we tr anslated t he observed J — H and H — K colors of 
stars ()Koornnee f 1983) to fluxes relative to K band and 
ignored the difference between K and Kg bands. We 
used a dust extinction model wi th a ratio of selective- 
to-to tal extinction (i?y) of 5.5 (jWeingartner &: Draind 
2001', their case "B"). Th is model is a good match to the 
Indebctou w et all (|2005f l extinction law study using the 
Spitzer IRAC bands. 

YSO candidates are traditionally selected from sources 
with excess flux at near- to far-infrared wavelengths com- 
pared to the stellar photosphere. Unfortunately, most 
(if not all) background galaxies fit the same description. 
In order to create a sample strongly enriched in YSOs, 
we have used several statisti cal criteria derived from the 
SWIRE data of ELAIS Nl (iSurace et al.l[200l to sepa- 
rate 'YSO candidates' from other objects with infrared 
excesses. The Spitzer Wide-area InfraRed Extragalac- 
tic (SWIRE) survey imaged six separate fields selected 
to be away from the Galactic plane and free from inter- 
stellar dust. The ELAIS Nl region, located towards the 
Galactic North Pole, presumably contains almost noth- 
ing other than non-YSO stars and background galaxies 
making it useful for understanding the distribution of 
infrared excess non-YSO sources. 

We used both IRAC and MIPS data together to iden- 
tify 'YSO candidates'. The criteria we used were empir- 
ically derived from SWIRE and the Serpens molecular 
cloud. Detailed discuss ion can be found in the deliv- 
ery documentation and iHarvev et al.l (|2007D . In brief, 
we constructed "probability" functions from three color- 
magnitude diagrams: [4.5] vs. [4.5] - [8.0], [24] vs. [4.5] - 
[8.0], and [24] vs. [8.0] - [24]. Based on where a source lies 
in each color-magnitude diagram, it is assigned a proba- 
bility of being a galaxy. The final probability is then the 
product of the three individual probabilities, and if it is 
less than a set value, the source is classified as a 'YSO 
candidate'. 

4. FINAL CATALOGS 

All catalogs discussed in the paper are available on the 
SSC website under the delivery 4 products'^. Separate 
catalogs of all 'YSO candidates' in Lupus are also avail- 
able from the same location. 

4.1. Main Lupus Catalogs 

^ i jhtt p : / /ssc . spit zer . caltech .edu/tools /starpetl 

^ http: / /ssc. spitzer. caltech.edu/legacy/c2dhIitory.html 



8 



In addition to the visual inspection performed during 
bandmerging, we made two further cuts on the data to 
obtain a high reliability catalog for each cloud region. 
First, we required all 24 /im detections to be at least 3a 
in the combined epochs (quality of detection "A" , "B" , or 
"C" ) . Asteroid contamination is also a problem in Lupus 
since the cloud is located near the plane of the ecliptic 
(latitude -14 to -21 degrees). To eliminate asteroids, we 
imposed the requirement that a source be detected at 
24 /im in both of the single epochs of observation and in 
the combined epochs data. Because our two epochs of 
observation were separated in time by 3 — 7 hours, the 
asteroids would have moved sufficiently that they would 
not be detected in both epochs at the same position. A 
broader statistical analysis of field as teroids in Spitzer 
Legacy d atasets will be presented in iStapelfeldt et al.l 
()in prep.h . Our catalogs contain 1790 sources in Lupus 
I, 1950 in Lupus III, and 770 in Lupus IV. In this paper, 
we only use these 'MIPS high-reliability' catalogs. 

4.2. c2d Processed 24 iim SWIRE Catalogs 

We have com pared our results t o the SWIRE data 
for ELAIS Nl (jSurace et al.ll2004D . Because this re- 
gion maps the Galactic north pole, the sources should be 
almost entirely non-YSO stars or background galaxies. 
Thus, this data set is valuable in distinguishing YSOs 
from galaxies. We processed a portion of the SWIRE 
BCDs for ELAIS Nl with the c2d data pipeline and cre- 
ated our own c2d-processed SWIRE catalog. This cata- 
log has only one epoch of observation for MIPS. However, 
ELAIS Nl is located at ~ -1-70 ecliptic latitude so aster- 
oid contamination is not a concern. 

The SWIRE data is more sensitive at all wavelengths 
than our c2d observations. To make comparisons be- 
tween SWIRE and c2d, we have created resampled 
SWIRE catalogs that simulate the sensitivity of our c2d 
observations. This resampling involves: (1) extincting 
the SWIRE sources consistent with positioning these 
sources behind each Lupus region; (2) randomly select- 
ing faint SWIRE detections to flag as non-detections in 
the resampled catalog in order to match sensitivity with 
c2d; (3) replacing the photometric errors with ones more 
characteristic of c2d observations; and (4) re-classifying 
the sources in the resampled catalog as described in § 13.41 

After resampling, we have three c2d-processed SWIRE 
catalogs, one each for Lupus I, III, and IV. To create 
'MIPS high-reliability' catalogs similar to above, we have 
also required a detection at 24 /im to be at least 3cr (qual- 
ity of detection "A","B", or "C"). Our catalogs contain 
1588 sources (Lupus I), 1646 (Lupus III), and 2061 (Lu- 
pus IV). Each one covers 5.3 deg^. We will use these 
resampled, 'MIPS high-reliability' c2d-processed SWIRE 
catalogs in this paper. 

4.3. Lupus Off- Cloud regions 

We observed two off-cloud regions, called OCl and 
0C2. For the analysis in this paper, we will only be 
using OCl. OCl is located at a higher Galactic lati- 
tude compared to 0C2. The Galactic latitude for OCl 
is closer to the latitude of the three cloud regions. Be- 
cause of this, the source counts and source classifications 
in OCl provide the best comparison with the on-cloud 
regions. We made the same 'MIPS high-reliability' cuts 
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Fig. 8. — Differential source counts at 24 microns for the 
three Lupus regions. Four different curves are displayed: The on- 
cloud source counts (black), Lupus off-cloud 1 (dotted), resampled 
SWIRE (gray), and 25 micron Wai nscoat galaxy model counts 
(dot-dashed) (Wains coat et ahl I1992f l . The Wainscoat model in 
each panel was computed for the Galactic latitude, longitude, dis- 
tance, and average Ay of the specified Lupus region using the 
values listed in table [2] 

to the data as used for the on-cloud catalogs. Our OCl 
catalog contains 157 sources. 

5. RESULTS 
5.1. 24 fim Source Counts 

Figure [5] shows the 24 /im differential source counts per 
square degree for Lupus I, III, and IV. For comparison, 
we also show the source counts per square degree for 
Lupus OCl, our c2d-processed SWIRE catalogs, and the 
predicted values from th e Wainscoat model at 25 microns 
(jWainscoat et al.lll992[ ). The Wainscoat model counts 
were produced separately for Lupus I, III, and IV using 
the Galactic longitude, latitude, distance, and median 
Ay present within the each region. These parameters are 
listed in Table [2l We did not apply a correction for the 
difference between the IRAS 25 fj,m counts output from 
the Wainscoat model and our plotted 24 fj,m differential 
source counts. The source counts for every curve were 
binned with a bin width of 0.1 mag., the same bin width 
produced by the Wainscoat model. 

In all areas, both on- and off-cloud, the counts peak 
around 1 mJy and drop off dramatically below this flux 
due to the completeness limit. Brighter than ^ 10 mJy, 
there are too few 24/im sources in OCl to make a reliable 
histogram. Therefore we have not plotted this curve past 
10 mJy. For fainter fluxes, the source counts for OCl are 
consistent with all three Lupus cloud regions. 

The predicted source counts from the Wainscoat model 
agree well with the actual on-cloud source counts for Lu- 
pus I, III, and IV. The only significant difference from the 
model is below ~ 3 mJy in Lupus I. Here, the on-cloud 
counts are higher than those predicted by the Wainscoat 
model, but are consistent with the resampled SWIRE 
counts, suggesting that these sources may be extragalac- 
tic. For Lupus HI and IV, the on-cloud counts are slightly 
higher than the SWIRE counts at all fluxes. Because Lu- 
pus HI and IV are consistent with the Wainscoat model, 
the higher counts are likely due to Galactic sources. This 
explanation is consistent with the relative Galactic lat- 
itudes of each cloud. Lupus I, at -1-17, should have a 
higher fraction of extragalactic sources compared to Lu- 
pus HI and IV which are at lower latitudes, and thus 
should have a higher fraction of Galactic sources. 

5.2. MIPS Color-Magnitude Diagrams 

We plotted color-magnitude diagrams for each Lupus 
region in Figures [9] and [TOl In each figure, the resampled 
SWIRE catalog for the specified Lupus region is shown 
with shaded contours. Since we did not create a resam- 
pled SWIRE catalog for Lupus OCl, we have shown the 
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TABLE 3 
Individual YSO Candidates 



Index Region Source Name RA (J2000) DEC (J2000) 



SSC ID 



Class 



1 


I 


IRAS 15356-3430 


15: 


:38: 


:48, 


,4 


-34: 


:40: 


:38, 


,2 


SSTc2d 


J153848, 


,4- 


-344038, 


,2 


I 


2 


I 


IRAS 15398-3359 


15: 


:43: 


:02, 


,3 


-34: 


:09: 


:07, 


,5 


SSTc2d 


J154302, 


,3- 


-340907, 


,5 


I 


3 


I 


2MASS 15450634 


15: 


:45: 


:06, 


3 


-34: 


:17: 


:38, 


1 


SSTc2d 


J154506, 


,3- 


-341738, 


1 


Flat 


4 


I 


2MASS 15450887 


15: 


:45: 


:08, 


,9 


-34: 


:17: 


:33, 


,7 


SSTc2d 


J154508, 


,9- 


-341733, 


,7 


II 


5 


III 


2MASS 16070854 


16: 


:07: 


:08, 


,6 


-39: 


14: 


:07, 


,7 


SSTc2d 


J160708, 


.6- 


-391407, 


,7 


Flat 


6 


III 


Sz 102 


16: 


:08: 


:29, 


,7 


-39: 


:03: 


:11, 


,0 


SSTc2d 


J160829, 


,7- 


■390311, 


,0 


I 


7 


III 


IRAS 16051-3820 


16: 


;08: 


:30, 


,7 


-38: 


:28: 


:26, 


,8 


SSTc2d 


J160830, 


.7- 


-382826, 


,8 


II 


8 


III 


IRAS 16059-3857 


16: 


:09: 


:18, 


1 


-39: 


:04: 


:53, 


,4 


SSTc2d 


J160918, 


.1- 


-390453, 


,4 


I 


9 


IV 


2MASS 16011549 


16: 


:01: 


:15, 


,6 


-41: 


:52: 


:35, 


,3 


SSTc2d 


J160115, 


.6- 


-415235, 


,3 


Flat 


10 


IV 


IRAS 15587-4125 


16: 


:02: 


:13, 


1 


-41: 


:33: 


:36, 


,3 


SSTc2d 


J160213, 


.1- 


-413336, 


,3 


I 


11 


IV 


IRAS 15589-4132 


16: 


:02: 


:21, 


,6 


-41: 


:40: 


:53, 


,7 


SSTc2d 


J160221, 


.6- 


-414053, 


,7 


I 


12 


IV 


2MASS 16023443 


16: 


:02: 


:34, 


,5 


-42: 


:11: 


:29, 


,7 


SSTc2d 


J160234, 


.5- 


-411129, 


,7 


Flat 



Note. — Names, positions, and classifications of individual YSO candidates discussed in ^15.51 The 
index numbers identify the sources in the figures. Classifications are based on the criteria of I Greene et al.l 



TABLE 4 

Fluxes for Selected YSO Candidates 







2MASS 






IRAC 






MIPS 




Index 


J 


H 




3.6 ^m 


4.5 


5.8 ^m 


8.UAtm 


24 ^m 


70 Atm 


160 ^m 




(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


1 


2.15 


3.84 


6.60 


3.95 


2.76 


9.58 


42.7 


114 


2920 


4230 


2 








1.97 


23.2 


42.6 


127 


961 


15300 


57200 


3 


0.955 


4.01 


9.61 


14.1 


15.4 


16.2 


19.1 


69.9 


204 




4 


21.1 


56.7 


87.4 


101 


104 


100. 


94.3 


137 


702 




5 


1.30 


5.73 


15.8 


32.9 


48.3 


60.6 


89.3 


196 


189 




6 


2.39 


3.54 


6.19 


14.4 


24.1 


34.0 


67.0 


347 


257 




7 


410. 


452. 


342. 










473 


2200 


2210 


8 








0.255 


1.00 


0.985 


0.549 


32.4 


2610 


8710 


9 


0.424 


1.67 


4.83 


8.36 


9.92 


8.98 


7.70 


75.9 


1220 


5390 


10 


2.94 


2.12 


1.97 




2.11 




3.80 


533 


276 




11 


1.63 


3.36 


4.66 


5.10 


4.87 


12.6 


42.1 


144 


1030 


2340 


12 


4.39 


7.45 


14.0 


31.1 




41.5 




118 


269 





Note. — The 2MASS and Spitzcr fluxes for the twelve selected sources in Table |3] The SEDs 
for these sources are plotted in Figure 1171 The statistical uncertainty in these fluxes is less than the 
absolute photometric error. As discussed in the delivery documentation IIEvans et al.ll2007h and in i) l3.2l 
the absolute photometric error is estimated to be 15% for 3.6 — 24 lira and 20% at 70 and 160 /^m. 



resampled Lupus I SWIRE catalog in the OCl panels. 
Overlaid on the contours are the data points from cat- 
alog for each Lupus region. We have plotted the data 
using one of four different symbols depending on our 
source classification: gray circles are stars, crosses are 
'YSO candidates', gray triangles denote sources outside 
of the observed IRAC area, and black triangles represent 
all other categories. Furthermore, in Figure [9] we have 
drawn a black circle around any source also detected at 
70 ^m. A few sources have been identified by number. 
These will be discussed in § 15.51 Note that source num- 
bers 2 and 8 are missing from Figure [9] because those 
sources do not have a detection at Ks but do appear in 
Figure [101 

Figure [9] shows the Kg vs. Kg — [24] plot for each cloud 
region and OCl. We have drawn do tted lines to separat e 
different YSO classes, as defined bv I Greene et alj ((1994). 
The four YSO classes are determined by the value of 
a — j°iog^^ 1 the slope of the linear least-squares fit to 
all data points from Kg ~ 24 /im. The four classes are as 
follows: Class I: a > 0.3; Flat spectrum: 0.3 > a > —0.3; 
Class II: -0.3 >a> -1.6; and Class III: a < -1.6. We 
have included sources with a = 0.3 into class I since these 



sources were left undefined bv lGreene eFall (|1994l ). We 
expect the SWIRE data should contain nothing but non- 
YSO stars and background galaxies (see i) l4.2p . Based 
on the shaded contours, we can group the sources into 
roughly one of three regions in this figure: stars are lo- 
cated at i^Tj- [24] « 0, galaxies are fainter than Kg ~ 13.5 
mag with Kg — [24] between 4-10, and YSOs have a 
Kg — [24] color excess and are brighter than Kg ~ 13.5 
mag. 

We have also plotted [24] vs. [24] — [70] for each region 
(Figure fTO]) . Only the brightest stars should be detected 
at lOfim and in fact we do see just one star, in Lupus III, 
with [24] - [70] « 0. Based on the SWIRE contours, we 
expect that sources with [24] > 6 mag are likely galaxies, 
while those brighter than this lim it are YSOs. We plot - 
ted four YSO model curves from IWhitnev et all (|2003f ) 
showing different stages of evolution. Each curve repre- 
sents a range of inclinations from pole-on (upper left) to 
edge-on (lower right). These models assume a distance 
of 140 pc, to correspond with nearby star-forming clouds 
such as Taurus. We shifted the 24 /Ltm magnitudes in 
each sub-panel to correspond to our assumed distance 
for that cloud region. We did not apply a shift to the 
off-cloud region. With only a few notable exceptions, 
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flat: Class I 




TABLE 5 
YSO Candidate Statistics 



10 12 2 4 6 a 10 12 
K,-{2'1] (mag) 

Fig. 9. — Ks vs. Ks — [24] color-magnitude plot for each Lupus 
region and the SWIRE region. The shaded contours are the co unts 
for the resampled SWIRE region, as described in the text in ii l4.2l 
The contour levels are 1,2,4,8,16, and 32 counts. The data points 
show the specified Lupus region: I, III, IV, or OCl. Different 
symbols are used based upon how a source is classified (^ I3.4| I. Gray 
circles are stars, crosses represent 'YSO candidates', gray triangles 
are sources only detected at JHKs and 24 fim wavelengths, and 
black triangles are all other classifications. Furthermore, a black 
circle surrounds any source also detected at 70 ^tm. The vertical 
dotted lines show the class boundaries as defined by Greene ct al. 
II1994I) . The numbers reference sources discussed in §[53] and listed 
in Tables [3] and H 



\ Lupus 1 


- 


\^ Lupus IV 
+ 10 4 \ ■■■^ 





Fig. 10. — [24] vs. [241 — [70] color-magnitude plots for each Lu- 
pus region. See Figure [9] for a description of the symbols, shaded 
contours, and numbe rs. Four protostar models are shown from 
IWhitnev etaLl ll2003l '). Except for Lupus OCl, we shifted the 
model tracks to correspond to our assumed distances for each Lu- 
pus region. 

our 'YSO candidates' are bluer than what is predicted 
by the models. 

5.3. Young Stellar Objects 

IComeroiil (j2007| ) has compiled the most up-to-date list- 
ing of all known and suspected YSOs in Lupus. Within 
our observed IRAC-j-MIPS areas, there are just 5 such 
objects in Lupus I, 55 in Lupus III, and 3 in Lupus 
IV. The number of 'YSO candidates' identified from our 
source classification mirrors this distribution; there are 
16 'YSO candidates' in Lupus I, 75 in Lupus III, and 
12 in Lupus IV. Compared to other c2d clouds such as 
Perseus, Ophiuchus, and Serpens, Lupus is relatively qui- 
escent; only of one the regions, Lupus III, has significant 
star formation activity. The number of 'YSO Candi- 
dates' per square degree is 12, 58, and 33, for Lupus 
I, III, and IV, respectively. Compared to these values, 





Lupus I 


Lupus III 


Lupus IV 


YSO candidates 


16 


75 


12 


YSOs per deg^ 


12 


58 


33 


Known YSOs (4")'' 


3 of 5 


27 of 55 


2 of 3 


Class 1 


5 


2 


1 


Flat Spectrum 


2 


5 


1 


Class II 


6 


45 


5 


Class III 


3 


23 


5 



Note. — Statistics on the number of 'YSO can- 
didates' identified in t he th ree Lupus regions based on 
IRAC and MIPS. See § l3.4l for a discussion of how these 
sources are defined. 

Number of m atches with kno wn and suspected YSOs 
(compiled from IComer6nir2007l) within 4 arcseconds. 

the number of ' YSO candidate s' per square degree is 
~ 300 in Serpens jHaryev et al.| [2006'). ~ 100 in Perseus 
(|j0rgensen et al. II2006D. a nd ^200 in the L1688 region 
of Ophiuchus (| Allen et al]|2007[ ). 

We have compared our 'YSO candidate' sample with 
the list of known and suspected YSOs to see how many 
of the sources we identify. Our results are summarized in 
Table [5] Overall, we find a match within 4" for 32 out of 
the 63 known and suspected YSOs. If we only consider 
previously known classical T-Tauri stars, we identify 26 
out of 40, or about two-thirds. Of the 31 "missed" YSOs, 
25 are not in our high-reliability catalogs while the re- 
maining six are not identified as 'YSO candidates'. In 
Table [5] we have also broken down our 'YSO candidates' 
by YSO class. For all three cloud regions, the majority 
of our objects are in the Class II and Class III categories. 
This result agrees with previous studies that have found 
a number of T-Tauri stars within the cloud, but no ev- 
idence of a large number of Class I and Flat spectrum 
stars fe.g. !KrauttCTlll991h . 

In Figures [TT] - [T3] we have plotted the positions of 
all our 'YSO candidates' on the 160 /im emission. Since 
we require both IRAC and MIPS data to classify 'YSO 
candidates', we have shown the extent of the IRAC ob- 
servations with white outlines. A few sources have been 
numbered; these will be discussed in the following sec- 
tion. 

5.4. Extinction Maps 

To create extinction maps for our observed regions, we 
started by excluding 2JVIASS sources with only an upper 
limit flux in any of the JHKs bands. For the remainder 
of the 2JVIASS sources, we calculated Ay values from the 
JHKs bands using the method outlined in ^ \'3A\ Our 
method gave us pinpoint measurements of the Ay along 
each source's line-of-sight. We converted these randomly 
distributed extinction values into uniform maps by creat- 
ing a 0.'4 pixel grid on top of the data. At each grid posi- 
tion, the extinction value in that cell equals the weighted 
average of the extinctions of individual points within a 2' 
radius. This average is weighted both by uncertainty in 
the individual Ay values and by distance from the center 
of the cell with a Gaussian function having a full-width 
half-maximum (FWHJVI) equal to 2'. Thus, our maps 
have 2' resolution with 5 pixels across the FWHIVI. We 
chose this cell size and resolution to have a well-sampled 
PSF and to ensure that our maps would have an extinc- 
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Fig. 11. — Lupus I 'YSO candidates' overlaid on the 160 fim 
image. The four symbols correspond to different YSO classes as 
noted in the legend and defined i n §15.21 A few sources have been 
numbered; these are discussed in § 15.51 Since our YSO classification 
procedures are based on IRAC and MIPS data, we have drawn a 
white outline to illustrate the area with both IRAC and MIPS 
observations. 

tion value for every pixel, even in the highest Ay regions. 
The median number of Ay values in each cell is roughly 
40 for Lupus I, and 80 for Lupus III and IV. The percent- 
age of cells with less than 10 Ay values is 1 — 2% for all 
three cloud regions. The extinction contours are shown 
in Figures [TH [151 HH overlaid on the 160 /xm emission 
(grayscale). We subtracted the median off-cloud Ay, 0.7 
mag., from each map. In all three maps, the extinc- 
tion contours start at 3 magnitudes of Ay and increase 
in steps of 3 mags. There is a strong correspondence 
between the 160 /J,m emission and the Ay contours sug- 
gesting the 160 /im emission traces the dust. A detailed 
comparison of Ay and 160 /im emission is beyond the 
scope of this paper. 

The mass in each Lupus region was estimated by con- 
verting the extinction value s in each cell t o colu mn den- 
sity using the equation from lBohlin et all (|1978f ). assum- 
ing Ry = Ay/E{B - = 3.1: 



NHi/Ay = 1.87 X 10^1 cm^ mag'^ 



(2) 



Note that this relationship was established for the dif- 
fuse ISM {Ry = 3.1) and not in dense molecular clouds 
such as we observed, so this equation may not be valid. 
This is an important concern since we have used the 
Weingartner & Draine Ry — 5.5 d ust model in thi s paper 
(see §133). The data compiled bv lKim fc Martini ()1996D 
are consistent with Nni/Ay being independent of Ry, 
albeit with large uncertainties. In this paper, we will use 
the value of Nni/Ay given in Equation [21 but note that 
this is possible source of error in our final mass estimates. 

With this equation, plus the molecular wei ght per hy - 
drogen molecule, fiH2 = 2.8 (calculated from lCoal2000f ). 
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Fig. 12. — Lupu s II I 'YSO candidates' overlaid on the 160 /im 
image. See Figure 1111 for a description of the symbols and their 
meanings. 

the mass of a hydrogen atom, and the distance to each 
Lupus region, we computed masses for each Lupus region 
from the extinction maps. In Table [2] we list the derived 
mass in regions with Ay > 3 for each cloud based on our 
extinction maps. 

It is difficult to compare our mass estimates with val- 
ues in the literature since the mapped areas are not al- 
ways comparable in size. However, two authors have 
observed comparable regions in ^^CO ( J = 1 ^ 0) and 
Ci^O (J = 1 ^ 0) (hereafter " CO and C^^O, respec- 
tively.) iTachihara et al.l (|1996[ ) observed Lu pus I and 
III in "CO with 8' resolution and lHara et"an ()1999[ ) ob- 
served Lupus I, III, and IV in C^^O with a resolution of 
2'6. Both papers computed cloud masses from emission 
> 3cr. With the exception of Lupus I, where these au- 
thors' maps extend about 1 degree further to the east, 
these > 3(T areas are similar to our Ay > 3 regions. For 
purposes of comparison to our results, we will scale these 
authors' mass estimates using our assumed distances to 
the Lupus clouds. This only involves changing their mass 
estimates for Lupus III. We have assumed a distance of 
200 pc to th is cloud whereas ITachihara et al.l ([1996) and 
iHara et all (1999.) assumed a distance of ~ 150 pc for all 
Lupus regions. 

The "CO mass in Lupus I is 1200 Mq which is ~ 3x 
larger than our extinction mass of 440 Mq. However, 
their mapped area is larger than ours and the 3cr ^^CO 
contours cover more area than our Ay > 3 maps. For 
Lupus III, the ^^CO mass is 530 Mq, very similar to 
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Fig. 13. — Lupu s IV 'YSO candidates' overlaid on the 160 fivci 
image. See Figure 1111 for a description of the symbols and their 
meanings. 



our extinction mass of 690 Mq . In Lupus I and IV, the 
C^^O masses agree well with our Ay values; the C^^O 
masses are 327 and 216 Mq, respectively, for Lupus I 
and IV compared to 440 and 250 Mq we derived from 
Ay. The only difference is in Lupus III, where our Ay 
mass, 690 A/q, is 3.5x that derived from C^^^O (ISTMq). 

Although some differences between our Ay mass esti- 
mates and those derived from ^'^CO and C^^O exist, there 
are several possible sources of error that must be consid- 
ered. Our extinction maps do not cover the exact same 
areas as the molecular cloud maps which can potentially 
lead to large differences in mass estimates, as is evident 
in Lupus I. Furthermore, the molecular masses were esti- 
mated assuming a constant excitation temperature of 13 
K, even though measurements of the excitation tempera- 
ture in individual dense cores in Lupus shown variations 
(jVilas-Boas. Mvers. fc Fulled [200 Q) . Similarly, our Ay 
masses assumed a constant grain composition through- 
out the clouds. This quantity ha s also been shown t o 
vary along different sightlines (e.g. ICardelli et aLlll989D . 
Furthermore, our Nhi /Ay ratio may also be incorrect. A 
detailed study is needed before making rigorous compar- 
isons between cloud masses derived using various molec- 
ular tracers and extinction measurements. 

5.5. Selected Sources 
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Fig. 14. — The magnitudes of visual extinction (contours) are 
overlaid on the 160 fim emission map for Lupus I. The contours 
range from 3 to 15 Ay, in steps of 3 Ay. The extinction map 
has 2 arcminute resolution. We have also shown the direction of 
increasing Galactic I and b coordinates. 

In this section we highlight a few sources we selected 
from Figures [51 and [TUl that are red relative to the general 
'YSO candidates' and hence hkely to be younger YSOs. 
From Figure [9] we selected all Flat and Class I sources 
bright enough in Ks to not overlap the shaded SWIRE 
contours denoting 'background galaxies'. Furthermore, 
we selected sources in Figure [TU] with [24] — [70] > 4 
and also not overlapping the shaded SWIRE contours. 
Our final sample consists of 12 objects. These sources 
are listed in Tables [3] and H] with positions, fluxes, index 
numbers, and YSO classification based on a. These in- 
dex numbers are also shown in the color-magnitude plots 
(Figures [9l-fT0|). in Figures [Tl] - [131 showing the location 
of all 'YSO candidate' sources, and in Figure [T71 where 
the SEDs for these objects are plotted. 

5.5.1. Lupus I Sources 

We selected four sources in Lupus I: one previously 
known and three objects without references in the liter- 
ature. The first object is IRAS 15356-3430. The relative 
brightness of this source in places it above the typical 
extragalactic region in Figure [51 This source is obviously 
extended in the 2MASS and IRAC bands. Furthermore, 
it is located in a relatively low extinction region to the 
south instead of near other known protostars. For these 
reasons, we believe it to be a background galaxy. 

IRAS 15398-3359 is a well- known protostar located 
in the core B228 (see, e.g. iHever fc GrahamI 119891 : 
I Shirley et all l2000l 120021 ). It has a molecular out- 
fiow which was discov ered in ^^CO (J = 2 ^ 1) by 
iTachihara et al.l (|1996[ ). There is some correspondence 
between the Herbig-Haro object HH 185 and the molecu- 
lar outflow: the blue-shifted velocity lobe from the molec- 
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Fig. 15. — The magnitudes of visual extinction (contours) are 
overlaid on the 160 fim emission map for Lupus III. The contours 
range from 3 to 18 Ay, in steps of 3 Ay. The extinction map 
has 2 arcminute resolution. We have also shown the direction of 
increasing Galactic I and b coordinates. 

ular outflow overlaps t he blue-shifted [S II] line from 
iHever fc GrahamI (| 19891 ). The [S II] Hue is expected to be 
excited by the shocks arising from the outflow jet inter- 
acting with the ambient medium. In the IRAC images, 
the protostar appears extended in the same direction as 
the molecular outflow. 

The SED of IRAS 15398-3359 is shown in Figure 
T7} The data poin ts at 450 and 850 fim come from 
Shirley et all (l2Ciol using their 40" (lower data points) 
and 120" (upper points) aperture values, while the 1.3 
mm point comes from Rcipurth ct al. (1993) with a 23" 
aperture. This source does not have a 2MASS detection 
so it does not appear in Figure [51 It is not classified as 
a 'YSO candidate' because it is saturated at IRACl and 
IRAC2 and therefore was rejected by our YSO classifica- 
tion procedure. 

To the east of protostar IRAS 15398-3359 is a cluster 
of objects with infrared excesses. We have plotted two 
sources from this region for which we could not find any 
references in literature. Both are near Sz 68, a known 
bright T-Tauri star. 2MASS 15450887-3417333 is - 50" 
away while 2MASS 15450634-3417378 is about - 80" 
away. 

5.5.2. Lupus III Sources 

Lupus III contains the most active star forming region 
in the Lupus cloud complex. HR 5999 dominates the 
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Fig. 16. — The magnitudes of visual extinction (contours) are 
overlaid on the 160 fira emission map for Lupus IV. The contours 
range from 3 to 27 Ay, in steps of 3 Ay. The extinction map 
has 2 arcminute resolution. We have also shown the direction of 
increasing Galactic I and b coordinates. 
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Fig . 17 . — SEDs of selected sources from the three cloud regions. 
See §[53] for a full discussion of these sources, and Tables [3l and 141 
for the positions and fluxes of the sources. We have included sub- 
millimeter and millimeter wavelength data from the literature for 
IRAS 15398-3359 and IRAS 16059-3857. For IRAS 153 98-3359, the 
data points at 450 and 850 /xm are from lShirlev et all (|200C1) with 
40" (lower point s) and 120" (uppe r) apertures. The 1.3 mm point 
comes from Rcip urth etal] l(l993|) with a 23" ap erture. The 1.2 
mm d ata point for IRAS 16059-3857 comes from FSchihara et al.l 
||2007| ). 

bright nebulous region in the south of our observed re- 
gion. This source is a well-known Herbig Ae/Be star. 
Because this source is saturated at 3.6, 4.5, and 24 fim 
wavelengths, we will not discuss it here. Three of our four 
sources are nearby, making HR 5999 useful as a location 
reference. 

2MASS 16070854-3914075 exhibits a fafrly flat SED 
indicating a strong IR excess over that of a stellar pho- 
tosphere. This source is located about 15' to the west of 
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HR 5999. There are several known YSOs in this region, 
but this source has not been noted in the hterature. 

Sz 102 is a weh-known T-Tauri star located just to 
the north of the nebulous region containing HR 5999. 
This sour c e is a lso known as Th 28 or Krautter's star. 
iKrautteJ (|1986f ) discovered an HH outflow from this 
source ( HH 228). There is some evidence for rotation of 
the jet (iCoffev et al. l [200l . The SED rises from 1.25 - 
24 fim, then drops at 70 /im. This source has also been 
observed by the c2d program with the Infrared Spec- 
trom eter (IRS) from '^5 — 35 /^m (|Kessler-Silacci et al.l 
l2006f l. The faintness of the object and the orientation 
of th e jet suggest that t he source may have an edge-on 
disk (|Hughes et al.lll99"^ . However, the c2d photometry 
and spectroscopy for this source do not clearly show the 
double-peaked str ucture that would b e expected for an 
edge-o n disk (e.g. IWood et al.ll2002D . IStapelfeldt et all 
(|2007f ) observed this source with HST and did not see 
any evidence for an edge-on disk. 

Our third source is IRAS 16051-3820. This object is 
outside of the observed IRAC area, so therefore it cannot 
be classified as a 'YSO candidate'. This source does have 
red Ks - [24] and [24] - [70] colors. In Figure UHl this 
is source #7 and is located between the Late Class I 
and Class II YSO models. In the 2MASS images, two 
objects, separated from each other by ~ 5" are clearly 
visible. These two objects are not resolved at 24 or 70/im. 
This source has not been previously identified as a 'YSO 
candidate' in the literature. 

Lastly, IRAS 16059-3857 is l ocated just to t h e east 
of the bright nebulous region. iNakaiinia "eTall (l2003h 
observed this source with JHKg imaging and discovered 
the region appears as a reflection nebula. They detected 
an elongated jet-like structure in Kg and a fan-shaped 
structure reminiscent of a cavity created by an outflow. 
The Herbig-Haro object HH 78 is just to the west. In 
the IRAC images, this source appears extended in the 
roughly the same direct i on as t he fan-shaped structure 
seen bv iNakaiima et all (|2003[ ). iTachihara eFall (|2007D 
found this source to be embedded in an H^^CO+ core and 
measured a 1.2 mm continuum PSF flux for it. Because 
there is no 2MASS data for this source, it does not appear 
in Figure However, it is reddest source in Figure [TUl 

5.5.3. Lupus IV Sources 

Compared to Lupus HI, Lupus IV is a more quiescent 
region with only a few T-Tauri stars, most of them pre- 
viously known. We focus on the four reddest objects 
in Figure [5] that are not within the 'background galaxy' 
region denoted by the SWIRE contours. 

2MASS 16011549-4152351 is located near the center of 
globule filament 17 (GF 17), at approximately the center 
of the mapped region. This source is consistently rising 
at all wavelengths longward of 8 /zm and appears to be 
a deeply embedded protostar in the filament. There is 
some nebulosity associated with the source at 3.6 and 
4.5 /im. 

IRAS 15587-4125 is the reddest source in Figure[Sl but 
has only a modest excess in Figure [TUl The SED shows 
a sharp rise between 8 and 24. This sou rce is th e known 
planetary nebula FN G336.9-f-08.3 (Ack er et al.] ir992). 

IRAS 15589-4132 has an SED that rises rapidly at 
wavelengths longward of 4.5 /im. Because this source 
is extended in all four IRAC bands and in the 2MASS 



JHKg images, we suspect this source to be a galaxy. 
This source has no identification in the literature. 

Our last source, 2MASS 16023443-4211294, is not iden- 
tified as a 'YSO candidate' because it falls outside of the 
observed area of IRAC2 (4.5/im) and IRAC4 {S^im). The 
only reference in the literature for this source is an x-ray 
source about 12" away. The source does not appear ex- 
tended in 2MASS, IRAC, or MIPS images. This source 
is located far from the globular filament, so it tempting 
to identify this object as a background galaxy. However, 
the SED for this object is fairly flat, which is consistent 
with other 'YSO candidates' instead of the rising SED 
found in our suspected galaxies. 

6. SUMMARY 

We presented maps of three regions within the Lupus 
molecular cloud complex at 24, 70, and 160 /Ltm observed 
with MIPS. We dis cussed the c2d data reduction pipehne 
(jEvans et al.ll2007l ) and the differences between the c2d 
pipeline and our catalogs used here. Our catalogs were 
bandmerged with 2MASS and SED fitting was performed 
to classify objects. After removing asteroids and select- 
ing only the highest reliability MIPS sources, our final 
catalogs contain 1790 sources for Lupus I, 1950 for Lupus 
HI, and 770 for Lupus IV. In addition, we have created 
three c2d processed catalogs from a subset of the ELAIS 
Nl SWIRE data. Each catalog was trimmed to have 
similar sensitivity and extinction as each of our cloud 
regions. 

We compared the 24 /im source counts in each region to 
our c2d processed SWIRE catalog, a catalog for a nearby 
off-cloud region, and the star counts derived from the 
Wainscoat model. We find that Lupus I, at the highest 
Galactic latitude, shows an excess of sources at fainter 
fiuxes over the Wainscoat model, while III and IV, at 
more modest latitudes, are consistent with Wainscoat 
over a broad range of fiuxes. The faint Lupus I counts 
are consistent with the c2d processed SWIRE catalog 
suggesting this region is more dominated by background 
galaxies. 

Usin g an empirical method developed bv lHarvev et al.l 
(|200l . we created a sample of objects that is strongly 
enriched in YSOs and minimizes contamination by back- 
ground galaxies. We identify 103 'YSO candidates' 
within the three Lupus regions surveyed. Of these, we 
find positional matches with 26 of the 40 confirmed clas- 
sical T-Ta uri stars in Lupus and 6 of the s uspected proto- 
stars f rom lLopez Marti et al.l (120051 ) and lComeron et all 
((200l . When we break down our 'YSO candidate' sam- 
ple by a, the best-fit slope to the flux from 2 — 24 /im, we 
find that almost all of the sources fall into the Class II 
and Class III categories. This suggests that the most of 
the protostars within Lupus are older and more evolved. 
There are a few of the younger Class I and Flat spectrum 
sources in each region. 

Our SED fitting procedure also computes the Ay for 
each object. We have used the 2MASS data to construct 
extinction maps of the three regions and compared them 
to the 160 /im emission. Our cloud masses derived from 
the extinction maps are largely comparable to previous 
studies, except for Lupus I, where previous authors have 
mapped a larger region than our observations, and in 
Lupus HI, where our extinction mass is 3.5 x larger than 
the C^^O mass from iHara et al.l (|1999f) . However, there 



15 



are several different factors that could alter these mass 
estimates such as assumptions about temperature or dust 
composition. A more detailed study is needed. 

Based on two color-magnitude diagrams we selected 
12 objects with very red Kg - [24] and [24] - [70] colors 
for individual mention. Three of these are previously 
known protostars: IRAS 15398-3359, Sz 102, and IRAS 
16059-3857. We suspect two other sources of being bright 
background galaxies, one object is a known planetary 
nebula, and one object is a double star. We are left with 
five previously unknown potential protostars. Further 
observations are needed to confirm whether these sources 
are bona- fide YSOs. 
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